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ABSTRACT 

IRAS04325+2402C is a low luminosity object located near a protostar in Taurus. We present new 
spatially-resolved mm observations, near-infrared spectroscopy, and Spitzer photometry that improve 
the constraints on the nature of this source. The object is clearly detected in our 1.3 mm interferometry 
map, allowing us to estimate the mass in a localized disk+envelope around it to be in the range of 
0.001 to 0.01 Mq. Thus IRAS04325C is unlikely to accrete significantly more mass. The near-infrared 
spectrum cannot be explained with an extincted photosphere alone, but is consistent with a 0.03- 
0.1 Mq central source plus moderate veiling, seen in scattered light, confirming the edge-on nature of 
the disk. Based on K-band flux and spectral slope we conclude that a central object mass > 0.1 Mq is 
unlikely. Our comparison of the full spectral energy distribution, including new Spitzer photometry, 
with radiative transfer models confirms the high inclination of the disk (> 80°), the very low mass of 
the central source, and the small amount of circumstellar material. IRAS04325C is one of the lowest 
mass objects with a resolved edge-on disk known to date, possibly a young brown dwarf, and a likely 
wide companion to a more massive star. With these combined properties, it represents a unique case 
to study the formation and early evolution of very low mass objects. 

Subject headings: stars: circumstellar matter, formation, low-mass, brown dwarfs - planetary systems 



1. INTRODUCTION 

The origin of brown dwarfs has been one of the most 
intensely debated subjects in cool stars research over the 
past decade. As of today, there is no consensus on the 
dominant mechanisms that determine the low-mass end 
of the stellar initial mass function. Recent modeling has 
provided predictions for observable properties in the very 
low mass (VLM) regime, which in principle would permit 
distinguishing between the proposed scenarios. In prac- 
tice, however, it has turned out to be challenging to place 
firm limits on the formation theory for VLM sources (se e 
reviews bv lWhitworth et al.ll2007tlLuhman et al.ll2007o[ K 

Of particular relevance for evaluating formation sce- 
narios are the global properties of disks and envelopes. 
The combination of high-resolution imaging and analy- 
sis of the spectral energy distribution (SED) in the in- 
frared/mm regime is generally acknowledged as the op- 
timum source of information for disk/envelope parame- 
ters of young stellar sources. As of today, constraints 
on radii and masses of brown dwarf disks are sparse: 
Recent deep mm observations provided first limits on 
disk ma sses (< 0.001 - 003M W ) and oute r radii (10- 
100 AU: IScholz et al J 12005 IKlein et al.lf2003h . Only one 
substellar disk has been resolved with high-resolution 
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imaging (|Luhman et al.ll2007af ). with a derived radius of 
20-40 AU. 

In this context, the young source IRAS04325+2402C 
(hereafter: IRAS04325C) in the Taurus star forming 
complex may be of particular interest. HST/NICMOS 
images of the Class I source IRAS04325+2402 
(|Kenvon fc Hartmannl [1995) . also known as L1535IRS, 
revealed IRA S04325C as a faint com panion, at a sepa- 
ration of 8'.'1 ()Hartmann et alj fl999) 6 . The object is re- 
solved 'into_a_d£u42kdobed_emission separated by a dark 
lane' (Hartma nn et al.lll999f ). resembling the images of 
edge-on T Tauri disks. Models including an edge-on disk 
and an envelope are able to reproduce the appearance of 
the images. From the near-infrared (NIR) photometry, 
the mass of this possible Class I object was estimated 
to be between 0.02 and 0.06 Mq, with an upper limit of 
0.25 M q. Based on statistical arguments. iDuchene et al.l 
(2004) find that the system is likely to be a physically 
bound binary (probability 98%). 

Thus, IRAS04325C is a potential benchmark object to 
study substellar formation. In this Letter, we present 
new millimeter interferometry and NIR spectroscopy as 
well as an analysis of the resolved SED for this source, 
including new Spitzer photometry, with the aim of im- 
proving the constraints on the properties of the sys- 
tem, particularly the mass of the central object and the 
disk+envelope. 

2. NEW CONSTRAINTS FROM RESOLVED 
OBSERVATIONS 

2.1. Millimeter interferometry 

Using the IRAM P lateau de Bure Interferometer 
(jGuilloteau et al.l [1992) , we obtained a 1.3 mm (226.50 
GHz) continuum map centered on IRAS04325C, ob- 

6 Note that the primary itself is likely to be a subarcsecond 
binary, see Hartmann et al. ( 199S|) 
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Fig. 1. — Interferometry map at 1.3 mm wavelength, obtained 
with the IRAM Plateau de Bure Interferometer. The map is cen- 
tered on our target IRAS04325C; the primary source is seen at 
(+2", -8"). The inlet in the lower left corner shows the beam size. 
The contours are 10% to 90% of the peak value in steps of 10%. 

served on 27 December 2007 (configuration C, rms 
0.3mJy/beam, HPBW 1.55" x 1.0"). The final map of 
128x128 square pixels with 0'.'53 pixel size was produced 
by Fourier transforming the calibrated visibilities using 
natural weighting using the Grenoble Software environ- 
ment GAG. 

The 1.3 mm image (Fig. [1]) shows two sources, which 
are identified with the primary (south) and companion 
(north ) object in the HST images from lHartmann et al.l 
( 1999). We estimate fluxes of 10±f mJy for both sources 
by integrating the pixel intensities. The 1.3 mm flux is 
attributed to thermal emission from cold dust. Since the 
peaks in the map are strongly centered at the positions of 
the infrared sources, the emission most likely originates 
from the disk and the localized envelope with a diameter 
< 300 AU. Assuming that the dust is optically thin, we 
can obtain an estimate of the total dust mass Md llst from 
the mm flux S v : 



-Mdust — 



(1) 



We assume a d istance D of 148 pc for the Taurus star 
forming region (|Loinard eraT][200l . a gas to dust ratio 
of 100, and a plausible range of values for dust temper- 
ature Td (10-25 K) and dust opacity k„ (l-3cm 2 g _1 at 
1.3 mm), as comm only suggested in the literature (see 
iScholz et al.ll2006l and references therein). Using these 
parameters, we obtain a total mass (dust and gas) of 1- 
10 Mj up and a most probable value of 3 Mj up for each of 
the two sources. This value is comparable to the high- 
est disk mass determined for a Class II brown dwarf in 
Taurus (jScholz et al.ll200ll . 

Single-dish measurements centered on IRAS04325 are 
available in the literature and can be used to put limits 
on the mass of a possible common envelope surrounding 
both the primary and the seconda ry. Published fluxes are 
~ 180 mJy at 850 ym (15" beam.lHogerheiide fc Sandell 
120001: I Andrews fc Williaml [2005|), 110 mJy at 1.3 mm 
CllMMotte fc Andrell200ll) . < 9.2 mJy at 3 mm (~ 11", 
lOhashi et al.l lT996). Using the same parameters as above 
and k v oc A _/3 with /3 ~ 1.5 (|Natta et al.ll2007ft yields 
corresponding masses of 10-30 Mj up of distributed dust 




Fig. 2.— IRTF/SpeX spectrum for IRAS04325C (solid line) in 
comparison with models (dotted lines). All spectra are scaled to 
the K-band flux. The models assume T e g = 3000 K and typical 
veiling. The three model spectra in the plot are calculated from 
three scattering models (MRN2, 12, 15 in lKenvon et al.N1993T) 

and gas in an area with diameter ~ 2000 AU around the 
primary IRAS source, including the companion. Fluxes 
integrated over larger areas (500 mJy for 60" diameter, 
iMotte fc And re 2001) clearly indicate the presence of 
large-scale mm background in that area. Given that the 
primary in the system is probably more massive than the 
companion (see Sect. [U, IRAS04325C is unlikely to be 
the dominant center of accretion in that area. Thus, the 
total mass reservoir for IRAS04325C is 1-10 Mj up in the 
localized disk/envelope plus presumably a small fraction 
of at most 30 Mj up in a possible common envelope. 

2.2. Near-infrared spectroscopy 

We obtained a NIR (1-2.5 urn) low-resolution (R ~ 
100) spectrum for IRAS04325C using SpeX at IRTF 
(|Ravner et al.l |2003[ ). Reduction, background correc- 
tion, and extraction was carried out using SpeXtools 
(jCushing et al.l I2004T) . We corrected for instrument re- 
sponse and telluric features using a spectrum of the AO 
star HD25175 observed immediately before the target. 
The final signal-to-noise ratio ranges from 15 in H- to 
30 in K-band. Our spectrum rises steeply from J- to ri- 
band and is lacking obvious photospheric features. The 
spectrum is shown in Fig. [2] (solid li ne) scaled to the 
published K-band flux (14. 0-14.6 mag, lHartmann et al.l 
Il999t IConnellev et al.ll2008h . In the following, we will 
use the shape of the spectrum and the K-band flux to 
derive a constraint on the mass of the central source. 

We compare the observed spectrum with models based 
on the phot ospheric templates STARdusty2000 and BD- 
dusty2000 (jAllard et al J 1200 If ). In a first attempt, we 
calculated a series of models assuming pure extinction, 
thus neglecting th e edge-on nature (NIR extinction law 
from iMathisll 199Clf ) . While it is possible to reproduce the 
slope of the spectrum using such models, the required 
combinations of T c g and Ay imply photospheric K-band 
fluxes that are more than one order of magnitude larger 
than the observed value. Thus, pure extinction cannot 
account for the NIR SED, confirming that the object is 
seen in scattered light through an edge-on disk. 

Thus, instead of extinction , we applied scattering mod- 
els from lKenyon et al. (1993), which have been confirmed 
to fit the colours of embedded young stars in Taurus. 
For a range of scattering properties and edge-on geom- 
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Fig. 3. — Spectral energy distribution for IRAS04325C in com- 
parison with Monte Carlo radiative transfer models. Triangles and 
diamonds are the observed fluxes; Spitzer datapoints at 3.6-8 fim 
are treated as upper limits. Solid line: model with ISM grains, 
dashed line: model with larger grains (see text), dash-dotted line: 
model with face-on geometry, dotted line: input photosphere. (All 
models include extinction of Ay = 10.) 

etrv. iKenvon et aH (|1993f ) obtain 2.0-4.5 mag difference 
between observed and photospheric K-band flux, in line 
with the properties of many known edge-on disks. This 
implies for our target Mk = 3.8 . . . 6.9 mag. Co mpari- 
son with the 1 Myr track by iBaraffe et alj (|1998l ) yields 
T cff = 2700.. .3200 K and M = 0.03 ... 0.16 M . 

Photospheric spectra for this range of temperatures ex- 
hibit deep H 2 absorption features, which are not seen 
in our spectrum. The most likely reason for their ab- 
sence is veil ing, as frequently ob served for young stel- 
lar sources ([Greene &: Ladal Il996| ). We use the typical 
veiling for Class-II sourc es of tk = 1-76 and rj = 0.97 
(jFolha fc Emerson! fl999h . If the K-band flux includes 
that amount of veiling, the upper limits for temperature 
and mass are reduced to T e g — 3000 K and M = 0.1 M Q . 

For both scattering and veiling, we fit the values for 
the photometric bands with a smooth function to ob- 
tain the wavelength dependence. Veiling is added first, 
assuming that all excess flux originates in regions close 
to the central object and is thus affected by scattering. 
Under these assumptions, we can fit the observed SED 
with T c ff ~ 3000 K, corresponding to a mass of ~ 0.1 M@ 
at 1.0 Myr. In Fig. [5] we show spectra f or three plau- 
sible scattering models (MRN2, 12, 15 in IKenvon et al.l 
1993), all with T e g = 3000 K and veiling, as described 
above. Stronger veiling, as u sually seen for Class I 
sources (|Doppmann et al1l2005f ). would further reduce 
the photospheric contribution to the K-band flux and 
allow for deeper H2O features, thus lower T e g and mass. 

Interpreting the NIR SED for IRAS04325C is chal- 
lenging and subject to large uncertainties. We have 
shown that under reasonable assumptions the observa- 
tional constraints are well-explained by veiled emission 
from a central object with 0.03-0.1 M & seen through an 
edge-on disk. Higher masses cannot be strictly ruled out, 
but would require unusual scattering characteristics. 

2.3. Spitzer photometry 

IRAS04325 has been observed with Spitzer 
IRAC/MIPS as part of the Taurus Spitzer Legacy 
Project (PI: D. Padgett). The spatial resolution of 
MIPS is not sufficient to properly resolve the system; 
the combined fluxes for the source, as given in the 



Taurus-1 source catalogue, are 1.9 Jy at 24 /im and 

4.7 ± 0.9 Jy at 70 fim. For comparison: The IRAS mid- 
infrared fluxes for this object are 2.1 Jy at 25 /im and 
12.86 Jy at 60 /im, again indicating that long- wavelength 
data with large apertures include significant extended 
background emission (see Sec. 12. ip . 

In the IRAC images, the object is clearly not point- 
like and shows two emission peaks coinciding with the 
positions of the two components, superimposed onto a 
bright nebulosity. We carried out resolved photometry 
at 3-8 /im by fitting the PSF in apertures with 3'.'6 ra- 
dius centered on the peaks of the two components. We 
estimate the fluxes of the companion to be 4.1, 3.5, 2.6, 

1.8 mJy at wavelengths 3.6, 4.5, 5.8, and 8.0 /im (uncer- 
tainty ±20%). Owing to the difficulties in disentangling 
the emission from the sources and from the underlying 
nebulosity, the true Spitzer fluxes of IRAS04325C are 
likely to be significantl y smaller. Comparing w ith the L- 
band flux published by lConnellev et all (2008) indicates 
that the background contamination may be as much as 
~ 70%. The fluxes of the primary source in the system 
are about 6, 7, 10, 11 times higher for IRAC channels 
1-4, respectively. 

3. MODELING THE SED OF IRAS04325C 

We compiled the SED for IRAS04325C from our new 
observations and the literature. In the N IR, we aim to re- 
produ ce the K- and L-band fluxes from lConnellev et al.l 
(2008). The Spitzer fluxes are considered to be upper 
limits, as outlined in Sec. 12.31 In Fig. [3] we compare 
the observed SED with three model SEDs from a Monte 
Carlo radiation transfer code, which has be en used suc- 
cessfu lly to model the SE Ds of T Tauri stars (IWood et al.l 
120021 ) and brown dwarfs (jScholz et al.ll2006h . The mod- 
els use NextGen m odel atmospheres for the photosphere 
(jAllard et al.ll200l[ ) and include dust destruction close to 
the central object. Accretion rates are assumed to be 
negligible for the heating of the disk. For more details 
on the model ingredients, see iScholz et al.l ((2006T ). In all 
three models, we fix th e disk radius at 30 AU, as given by 
lHartmann et a l. (1999) based on the analysis of the HST 
images, the extinction at Ay = 10, and the parameters 
of the central object at values consistent with our spec- 
troscopy (T off = 2600 K, M = 0.07 M Q ). In addition, 
the model contains an infalling envelope with a bipolar 
cavity. 

The models shown with solid and dashed line assume 
an inclination of i = 80° and two different grain size dis- 
tributions: small ISM-type grains (solid line) and larger 
grains (dashed line) , mimicking the effect of dust coagu- 
lation. In the latter case, the grain size distribution is a 
power law with an exponential decay for particles with 
sizes abov e 50 /xm and a form al maximum grain size of 
1 mm, see IWood et al.l (|2002f ). Both models are able to 
provide a good match to the available datapoints. To 
fit the mm datapoint we need a total circumstellar mass 
of 2.2 x 10" 3 M Q for ISM grains and 4.6 x 10~ 4 M Q for 
the larger grains. In the second case, the value for the 
disk/envelope mass is lower than the constraint given in 
Sec. 12.11 due to the higher dust opacity. Note that the 
two models are highly discrepant in their prediction for 
the fluxes at 5 — 20 fim due to the different scattering 
albedo of the two grain types. 

For illustration purposes, we also provide the SED for a 
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face-on geometry in dash-dotted lines. Clearly, the model 
is not in agreement with the observations. No combina- 
tion of T c ff and Ay is able to fit both the observed fluxes 
and the near-infrared spectral slope, as discussed in Sec. 
I2~2l 

In summary, the modeling clearly confirms the plausi- 
bility of the results obtained in earlier sections. The ob- 
served SED can be reproduced assuming that the central 
source is a very low mass star or brown dwarf, the incli- 
nation is high (i.e. the object is seen mostly in scattered 
light), and the circumstellar mass is negligible compared 
with the object mass. The models also demonstrate that 
more complete wavelength coverage is needed for further 
analysis of the properties of the system. 

4. THE NATURE OF IRAS04325C 

We present new observational constraints on the nature 
of the potential young brown dwarf IRAS04325C, a 8'.'1 
companion to a well-k nown IRAS source discovered by 
lHartmann et all (| 1999( 1 in HST/NICMOS images. From 
our NIR spectrum for IRAS04325C, we estimate that the 
central source likely has a mass between 0.03 and 0.1 Mq. 
The available constraints on the infrared SED are incon- 
sistent with a reddened photospheric spectrum, but they 
can be understood by assuming that the object is seen in 
scattered light. This is independent confirmation for the 
edge-on geometry of the disk, as seen in the HST images. 
IRAS04325C is only the second case of an edge-on sys- 
tem in the mass range < 0.1 M© that has been resolved 
with high-resolution imaging. 

We detect the source with 1.3 mm interferometry at 
a flux level of lOmJy, corresponding to a mass of 1- 
lOMjup in a localized disk/envelope with R < 300 AU. 
Since the accretion reservoir in a larger-scale envelope is 
probably insignificant, the central source is unlikely to 
gain more mass than 0.01 M© due to further accretion. 
Thus, IRAS04325C will probably become either a very 
low mass star close to the substellar boundary or a brown 
dwarf. 

The evolutionary stage of IRAS04325C is not well- 
determined based on the current data. The primary 
component of the system has been classified as a Class 
I source bas ed on the mid-infrare d spectral slope (a(2 — 
25) = 0.79, iDucheneet all I2004D and the NIR colours 
(lLada fc Adamslll992fl . It is also assumed to drive a CO 
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outflow (e.g. lHogerheiide et al.lll9 98). another character- 
istic for a protostellar source. The companion is clearly 
embedded in the same cloud core and likely physically 
bound to the primary, thus the two sources can be as- 
sumed to be coeval. This would make IRAS04325C the 
lowest mass Class I source known to date. The classi- 
fication of the primary source, however, is partly based 
on unresolved data, including light from the companion 
and the surrounding nebulosity. Furthermore, the disk 
of the primary is p robably seen close to edge-on (i ~ 80°, 
IFurlan et all [2008). which would suppress the NIR flux 
and might lead to its m isidentification as Class I source 
(see I White et al J [2007( 1. We estimate a total mass of a 
few Jupiter masses in a localised disk+envelope for the 
primary source, which seems to be too low for a Class 
I object (| Andrews fc William si [2001 . Thus, the Class I 
nature of the source is questionable. 

The central source of IRAS04325C has probably sig- 
nificantly less mass than the primary component in the 
system. In all near /mid-infrared-bands, IRAS04325C is 
about an order of magnitude fainter than the primary 
source, indicating a mass ratio of ~ 0.15. The luminos- 
ity of the primary has been estimated to be 0.7-0.9L© 
(|Kenvon et all 19931 ; IFurlan et ai1l2008f) . implying a mass 
of ~ 0.7 Mq. IRAS04325C also has a much smaller disk 
than the primary object, as seen in the HST images. 

With its combined properties - very low mass central 
object, edge-on disk, wide companion to a more massive 
star - IRAS04325C is a prime target for follow-up studies 
and presents us with a unique opportunity to study the 
early evolutionary stages in a critical mass range. 
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